Total internal reflection and evanescent gain 
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Total internal reflection occurs for large angles of incidence, when light is incident from a high- 
refractive-index medium onto a low-index medium. We consider the situation where the low-index 
medium is active. By invoking causality in its most fundamental form, we argue that evanescent 
gain may or may not appear, depending on the analytic and global properties of the permittivity 
function. For conventional, weak gain media, we show that there is an absolute instability associated 
with infinite transverse dimensions. This instability can be ignored or eliminated in certain cases, 
for which evanescent gain prevails. 

PACS numbers: 42.25.Gy, 42.25.Bs, 42.70.Hj 



I. INTRODUCTION 

When light is incident from a high-refractive-indcx 
medium onto a low-index medium, it undergoes total in- 
ternal reflection provided the angle of incidence is larger 
than a certain critical angle. Total internal reflection is a 
fundamental physical phenomenon with several famous 
applications; in particular modern telecommunications 
rely on optical fibers based on this phenomenon. 

Since the tangential electric and magnetic fields must 
be continuous at the interface, there must be nonzero 
fields in the low-index medium, even though the inci- 
dent wave is totally reflected. For lossless/gainless me- 
dia, these evanescent fields decrease exponentially away 
from the interface. The presence of evanescent fields in 
the low-index medium suggests that the reflected wave 
will sense any perturbation induced there. In particular, 
if the low-index medium has gain, the reflection response 
will change compared to the lossless/gainless case. The 
problem of determining the correct electromagnetic re- 
sponse in the case of an active low-index medium is far 
from trivial, and has been discussed for 40 years without 
reaching consensus 

GHHI. A key issue is whether the 
reflectivity may exceed unity (i.e., evanescent gain ex- 
ists) when the active medium fills the entire half-space. 
Exp eriments have indicated that evanescent gain exists 
[13 - 1151 ] . However, it has been argued that the ampli- 
fied reflection may be due to backreflection from e.g., the 
boundaries of the active medium [ll[ . 

When the active medium has a finite thickness, it is 
well known that the overall reflection from the slab may 
exceed unity. This situation is fairly simple, as there is no 
need to determine the sign of the longitudinal wavenum- 
ber in solving Maxwell's equations for this case; the two 
waves (with opposite signs) are present simultaneously. 

Since there are no gain media with infinite thickness, 
why examine this case? The answer becomes clear if 



we formulate a similar question in terms of the refrac- 
tive index: Why define the refractive index as a separate 
parameter, when the electromagnetic field in any realis- 
tic, bounded structure can be expressed in terms of the 
permeability and permittivity? While the refractive in- 
dex or longitudinal wavenumber are not needed to ob- 
tain the formal solution to Maxwell's equation in a finite 
slab, it is still useful since it immediately provides in- 
formation about the involved physics. For example, it 
predicts whether the medium refracts positively or neg- 
atively [16j . Also, assuming darkness for time t < 0, the 
solution to Maxwell's equations for a semi-infinite gain 
medium equals that of a finite slab for times t less than 
d/c, where d is the slab thickness and c is the vacuum ve- 
locity of light. Hence, understanding semi-infinite media 
helps explaining transient phenomena. 

We will now summarize the existing controversy. As- 
suming well defined frequency-domain fields, Maxwell's 
equations can be solved in the frequency-domain, using 
the sign convention exp(— iut). With respect to Fig. Q] 
we define the transverse wavenumber (spatial frequency 
of the source) k x . For simplicity we assume both media 
to be nonmagnetic. Let t\ and t 2 be the relative permit- 
tivities of the high-index medium to the left and the low- 
index medium to the right, respectively. For plane waves, 
Maxwell's equations require the longitudinal wavenum- 
bers in the high-index and low-index media to be 



hz = ±y / e 1 uj 2 /c 2 - fcf, 
k 2z = ±\Je 2 Lu 2 /c 2 - k 2 x . 



(la) 
(lb) 
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At some observation frequency uj = we assume 
k 2 < Reeiw 2 /c 2 while k 2 > Ree 2 w 2 /c 2 . Since the high- 
index medium is passive, we may readily determine the 
correct sign of the square root in Eq. (fTa|) . For the low- 
index medium, we assume Im e 2 < and | Im eg| <SC 1 
(i.e., small gain). The correct sign for the square root in 
Eq. (|lb[) is far from obvious: Either Imk 2z > and 
Rek 2z < 0, or lmk 2z < and Rek 2z > 0, see Fig. 
[2] None of these solutions are appealing: The first re- 
quires the phase velocity and Poynting vector to point 
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FIG. 1. A wave is incident from a high-index medium to a low 
index medium with gain. The source produces a single, spa- 
tial frequency k x . The electromagnetic boundary conditions 
require preservation of the wavenumber k x parallel to the in- 
terface. The longitudinal wavenumbers are denoted ki z and 
fez- Note that since the excitation is assumed to be causal, 
it contains a band of frequencies, and therefore also a band 
of fci z 's and fez's. 



towards the boundary. Since there are no sources at 
z = oo, one may argue that this scenario cannot be 
true [Tl|. The second solution requires that the fields 
increase exponentially away from the boundary. Also, in 
the limit of zero gain the fields will increase exponen- 
tially as exp(zy/k^. — Ree2Wj ! /c 2 ) (see Fig. [2|, while in 
the limit of zero loss, the fi elds decrease exponentially as 
cxp(— Zy/hZ — Ree2wf/c 2 ). Such a discontinuity seems 
unphysical [9]. 



In this work we will first go back to fundamental 
electromagnetics, to ensure that we use the principle of 
causality in its most primitive form: No signal can prop- 
agate faster than the vacuum velocity of light. After the 
general analysis in Section UTJ we consider conventional, 
weak gain media in Section IIIII and show that they pro- 
vide evanescent gain. In Section IIVI we present an ex- 
ample that demonstrates that not all gain media give 
evanescent gain; this depends on the medium's global 
dispersion behavior. 



Im fco 
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FIG. 2. The two possible solutions for the wavenumber fe z 
for monochromatic analysis and a gainy medium. The arrows 
indicate the two possible wavenumbers in the complex plane, 
as the gain tends to zero. For a lossy medium, we always have 
a solution that tends to the upper alternative +i\/k% — uj\jc 2 
in the limit of zero loss. For simplicity we have taken Re £2 = 1 
here. 



II. LAPLACE TRANSFORM 
FREQUENCY-DOMAIN ANALYSIS 



Going back to fundamental electromagnetics, we note 
that Maxwell's equations, combined with appropriate 
causal constitutive relations, contain everything neces- 
sary to obtain a unique solution. To determine the cor- 
rect solution, we must be certain that we consider the 
real, physical situation. The real physical fields are the 
ones in the time-domain. By requiring the fields to be 
zero for t < (see Appendix IB]) . we obtain the causal so- 
lution to Maxwell's equations. The complex frequency- 
domain fields are usually found from the time-domain 
fields by a Fourier transform. However, when there is 
gain in the system, using the Fourier transform can be 
perilous, since the field may increase with time. At first 
sight, any instability seems to be convective in our case. 
This is however not true: A causal excitation involves 
an infinite band of frequencies. For a single spatial fre- 
quency k x this means that modes with a wide range of 
incident angles are involved; in fact even the mode with 
k2 Z = may be excited. This "side wave" gets ampli- 
fied and leads to infinite fields at the boundary. This 
instability is somewhat artificial, since its existence is 
dependent on infiniteness in the transverse direction; we 
will argue below how it can be ignored in certain situa- 
tions. Nevertheless, within a linear medium framework, 
Fourier transforms do not necessarily exist. Therefore, 
as in electronics and control engineering, we generalize 
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the analysis by using the Laplace transform, 



E(oj) 



£(t) exjp(iu>t)dt. 



(2) 



In Eq. ([2]) a sufficiently large value of Im uj will quench 
an exponential increase in the time-domain electric field 
£(t), such that the integral converges. (Note that uj is 
complex in general, equal to is, where s is the conven- 
tional Laplace variable.) The inverse transform is given 

by 



£m = rJ_ 



+ 00+27 



E{oj) exp(— iujt)duo. 



(3) 



00+27 



The integral is taken along the line uj = 17, for a suf- 
ficiently large, real parameter 7, above all non-analytic 
points of E(uj) in the complex w-plane. An important 
observation is the following: The frequency-domain field 
E(uj) only has physical meaning through the transforms 
©-([3]). Thus, if the field is to be interpreted for all 
real frequencies, it must be analytic in the upper half- 
plane Imuj > 0. However, as is shown below, if the non- 
analytic points are located in the upper half-plane, but 
close to the real axis and far away from the excitation 
frequency, we can still attribute a physical interpretation 
to the frequency-domain expressions. 

In order to derive the Fresnel equations and deter- 
mine the sign of k 2z , it is tempting to start with the 
response from a slab of finite thickness d, and then take 
the limit d — > 00. For finite d the solution to Maxwell's 
equations is independent of the sign of k 2z in the slab 
[13, 13- However, for an active slab the multiple reflec- 
tions may diverge, especially for a large d. Thus, for real 
fre quen cies, the limit d — ^ 00 is not necessarily meaning- 
ful mini, a way around this, is to evaluate the fields 

for sufficiently large Imw, where the frequency-domain 
fields exist. There, an exponential increase is quenched 
by the exponential factor exp(— Ima; t). As a result, we 
can take the limit d — > 00 [16j ■ For TE polarization, the 
Fresnel reflection coefficient p and the transmission co- 
efficient r (including the propagation factor exp(ik 2z z)) 
become [3, Ojj 



P = 



-<2z 



kiz + hz 
2k lz 

kiz + k 2z 



exp(ifc 2z -z), 



(4a) 
(4b) 



provided the sign of k 2z is determined such that k 2z — > 
+uj/c as Imw — > 00, and k 2z is an analytic function of uj. 
Indeed, even though Eqs. ((4]) have been derived for large 
Imw, we can extend their valid region as follows: The re- 
flected and transmitted frequency-domain fields are given 
by Eqs. (g]) multiplied by the Laplace-transformed inci- 
dent field. The associated, physical, time-domain fields 
are obtained by the inverse transform ([3]). Now, by an- 
alytic continuation, we can reduce 7 until we reach a 
non-analytic point of Eqs. Q, without altering £(t). If 



the expressions ([4]) are analytic in the entire, upper half- 
plane, we can set 7 = and interpret p and r for real 
frequencies. On the other hand, if there are non- analytic 
points in the upper half-plane, the time-domain fields di- 
verge. In that case, real frequencies are not physically 
meaningful in general. 



III. WEAK GAIN MEDIA 

To find the actual reflection and transmission response, 
we first consider conventional weak gain media, with the 
following assumptions or properties: 

1. The permittivity e 2 (w) obeys the Kramers-Kronig 
relations. 

2. The gain and dispersion is small, so that the per- 
mittivity can be written 



e 2 (uj) = e 2 + Ae 2 (uj), |Ae 2 (w)| < e 2 



(5) 



for real uj. Here e 2 is required to be a positive con- 
stant. In the following we take e 2 = 1; the analysis 
can easily be generalized to the case with another 
e 2 . (In the latter case, e 2 is only constant in a wide 
frequency band including the band where Ae 2 (w) 
is nonzero; for very high frequencies it necessarily 
tends to 1.) 

3. The medium is gainy at the observation frequency 
u>i and the critical frequency k x c. 

4. Let Ae max = max^ |Ae 2 (u;)|. In a bandwidth 
AcmaxfcrC around the critical frequency k x c, the 
permittivity e 2 (w) varies slowly: 



de 2 
duj 



< — — for \w — k a 



c < 



(6) 



Properties 2 and 4 essentially mean that the gain is weak 
and the dispersion is small. 
We now solve the equation 



£ 2(wW 



(7) 



to determine whether k 2z has branch points in the upper 
half- plane of the complex oj plane. Since e 2 (w) satisfies 
the Kramers-Kronig relations, it is analytic in the upper 
half-plane. The maximum modulus principle of complex 
analysis [l9| therefore ensures that property 2 is valid also 
in the upper half-plane, not only at the real frequency 
axis. Substituting e 2 (oj) = 1 + Ae 2 (u;) into Eq. (JT)) we 
find 



CO 



±k x c 1 



Ae 2 (w) 



(8) 



in the upper half-plane, since |Ae 2 (w)| <C 1. Thus, every 
solution to the dispersion relation in the upper half-plane 
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FIG. 3. The complex u;-plane. For conventional, weak gain 
media, there are branch points right above oj = ±(k x c + S') ~ 
±k x c. The branch cuts can be chosen arbitrarily; however, the 
shown, vertical cuts minimize the integral around the part of 
the branch cuts in the upper half-plane. 



is located within a distance (Ae max /2)fc :c c from the crit- 
ical frequencies ik x c. 

We therefore examine the region around k x c in more 
detail. If there were two solutions uj a and uji, to the dis- 
persion relation, then Eq. would predict that 



Ae 2 (w a ) - Ae 2 (w 6 ) = ^( w 6 - w a ). (9) 

By property 4 this is impossible unless Ub = uj a - Thus 
there is a unique solution to Eq. ([7]) in the first quadrant, 
located in the vicinity of k x c: 



uj = k x c + 5' + id, 



(10) 



where 



S' = -^^l kx c aad^=-5^Mfc,c. (11) 

In addition there is a mirrored solution in the second 
quadrant, located at w = — k x c — 5' + i<5. Note that 



e , A6 max 

o < — - — k x c. 



(12) 



In the expression for fc 2z and the Fresnel coefficients 
(|4|, these solutions appear as branch points. Hence, 
when evaluating the physical time-domain fields by the 
inverse Laplace transform, we must integrate above the 
associated branch cuts, from —00+17 to +00+17, see Fig- 
El By path deformation this path is the same as the path 
from —00 to 00 plus the paths around the branch cuts 
in the upper half-plane (Fig. [3]). Thus we may use the 
inverse Fourier transform to determine the time-domain 
fields, but only if we add the integrals around the branch 
cuts. Due to the exponential factor exp(— iuji), the inte- 
grals around the branch cuts diverge and dominate after 
some time. 

The divergence of the time-domain fields can be ex- 
plained as follows. Any causal excitation involves an infi- 
nite frequency band. For example, the Laplace transform 
of a unit-step- function modulated cosine, u(t) cos(ujit), 
is iu)/(ui 2 — loI). Thus, it is nonzero for all finite w^O. 
One of these frequencies is the branch-point frequency 



for which fc 2z = 0, that is, uj ~ k x c + id. This frequency 
is complex; the imaginary part 5 means that the associ- 
ated eigenmode is a growing wave with envelope exp(<5t). 
Physically, a wave with fc 2z = propagates along the 
boundary. Because the medium is gainy, this side wave 
picks up gain on its way. Consider a fixed observation 
point, e.g. the point z — + and x — 0. Since the 
medium and the excitation are unbounded in the trans- 
verse ^-direction, there are side waves that start arbitrar- 
ily far away from the observation point. Thus the field at 
the observation point diverges. As the field in medium 
2 becomes infinite, the field in medium 1 is infinite as 
well. Since the field at a fixed point in space diverges 
and the instability is not a result of amplified, multiple 
reflections, the instability for the system in Fig . Q]can be 
classified as an absolute instability [l?], l2l| . 

This instability could be eliminated (or converted into 
a convective instability) by limiting the extent of the gain 
medium in the transverse direction with an absorbing 
boundary. Alternatively, the incident wave itself could 
be limited in the x-direction, leading to an infinite spec- 
trum of k x modes (see Appendix A and Ref. [H). Rather 
than imposing such remedies, we will simply calculate the 
time-domain fields by an inverse Laplace transform above 
the branch cuts. If the excitation frequency u>± is suffi- 
ciently remote from the branch points, the side wave with 
&2z = is only excited very weakly, and can be neglected 
up to a certain time. The condition that the excitation 
frequency is remote from k x c means that the incident 
angle is not close to the critical angle. This condition is 
imperative in order to distinguish between the reflected 
wave, with an angle of reflection equal to the angle of 
incidence, and the wave associated with the growing side 
wave, with "reflection" (or propagation) angle equal to 
the critical angle. 

The reflected time-domain field for the excita- 
tion u{t) exp(ik x x — iuiit), with Laplace transform 
exp(ik x x) / (iuii — iuj), is given by 



£ p (x,t) 



1 r+°°+ l i _ f t2z CX p(ik x x — iuot) 



2tt 



00+27 



ku + k 2z 



do;, 



ZWl 



(13) 

at z = 0. The integral (IT31 can be evaluated by a general- 
ized version of the residue theorem, in which we find the 
contour integral around all poles and branch cuts of the 
integrand in half-plane Imw < 7. Provided lo\ is suffi- 
ciently remote from any resonances of the two media, the 
transients due to all poles and branch cuts for Imcj < 
can be ignored. Alternatively, for times larger than the 
maximum inverse bandwidth of the resonances, the 
transients will have died out. Then the reflected field for 
x = is given by 



S p (0,t) = 



ku 



: <2z 



■expH^+fbctO,*), t>T~\ 



(14) 

where the wavenumbers k\ z and fc 2z have been evaluated 
at the frequency uj\. The term £^ c (0,t) is the integral 
(|13l) around the two branch cuts above lo = ±k x c. This 



FIG. 4. The function fclz(w) — ^i(uj)uj 2 j c 2 — k 2 for a typical 
gain medium, plotted in the complex fclz-plane. To identify 
&2z, we require it to be +uu/c at w = oo, continuous as uj 
decreases towards zero, except at the branch cut at uj = k x c 
where it changes sign. 



integral is bounded by 

|£bc(o,t)| < 



const 
k x c — U>l 



exp(Ji). 



(15) 



Here, the constant depends on the specifics of the active 
medium (see Appendix C). In other words, for < 
t ^ an d provided u\ is not too close to k x c, we can 
ignore £bc(0, t). Then the reflected field is well described 
by the first term in Eq. (fl"4|) . 

We can now answer the question about the existence 
of evanescent gain. To obtain Eq. (O, we have only 
considered two branch cuts in the upper half-plane; 
these are the necessary branch cuts due to the zeros of 
e 2 (ui)uj 2 /c 2 — k 2 . We must ensure that the integrand 
in Eq. (|13|) is analytic everywhere else in the upper 
half-plane. That is, the sign of k 2z must be determined 
such that k 2z is analytic everywhere, except at the two 
branch cuts in the upper half-plane. Since k 2z — >■ +uj/c 
as uj — > +00, we can determine the sign by decreasing u> 
from +00 to u>i, ensuring that k 2z is continuous every- 
where except at uj = k x c where it changes sign. From Fig. 
S] we find that Im k 2z > at the observation frequency 
u)\. Hence, for weak conventional gain media, provided 
the "reflected" field from the side wave can be ignored, 
evanescent gain is possible. This result is consistent with 
0, H| , and the time-domain simulations in [Toj | where the 
dispersion of the medium is discarded. 

In Fig. [5] we plot the reflected and transmitted elec- 
tric field for a weak Lorentzian medium, after the tran- 
sients have died out, and before the side wave dominates. 
The reflected field was computed by Eq. (fl"3")) . including 
the propagation factor exp(— ik\ z z) in the integral. The 
transmitted field was computed with the same equation, 
but with r instead of pexp(— ik\ z z) in the integral (see 
Eq. (01). For z > we clearly see an evanescent decay- 
ing field, while the reflected field for z < is larger than 
unity. 

It is interesting to examine the situation when we ap- 
proach the critical angle associated with the frequency 
ui\. If we insist on using only the first term of Eq. (fl"4)) 
in this case, a simple calculation shows that the power 
reflectance would have been bounded by (y/2 + 1) 2 ~ 5.83 
at the critical angle. Also, the wavenumber k 2z and the 
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FIG. 5. The reflected electric field (solid line, z < 0) and 
transmitted electric field (dashed line, z > 0) for a plane 
wave incident to a weak Lorentzian gain medium (|C1|I . The 
parameters used: F — 0.01, V — O.luo, k x c — 2uj a , wj = uj a , 
7 = 0.001, and ej = 4.7. The field is plotted for x = 
and uJot = 10 3 , and normalized to the incident field. The 
amplitude of the reflected field is 1.01. Note that the field 
is discontinuous at z = because the incident wave is not 
included. 



reflected field would be discontinuous as we pass the crit- 
ical angle. This is clearly a paradox, as the branch cuts 
were chosen arbitrarily. The dilemma is resolved by not- 
ing that the entire Eq. (| X4[) must be used in this domain; 
both terms naturally coexist and cannot be separated. 
As we approach the critical angle, £bc(0, t) becomes com- 
parable to or larger than the first term in Eq. (TT4"|) . for all 
times. A different choice of branch cuts will alter each of 
these contributions, but the sum remains the same. For 
finite transverse dimension, the side wave's contribution 
to the "reflected" field does not necessarily diverge any 
more; however, the intensity of the reflected field can be 
arbitrarily large as the dimension is increased, or if the 
reflections from the transverse end facets are large. 



IV. GENERAL GAIN MEDIA 

More sophisticated gain media can be constructed, at 
least in principle, that behave differently compared to the 
conventional weak gain media. We will here show that we 
can obtain a near-imaginary k 2z with negative imaginary 
part at an observation frequency oji < k x c. Consider the 
permittivity 



£2(w) 



(uj - N)(uj + N*) 
(uj -P)(uj + P*) 



uj\ 



(16) 



where the complex numbers N and P are located in the 
lower half-plane, and uj 2 is a real constant. The longi- 
tudinal wavenumber satisfies k 2z (uj) — e 2 (uj)uj 2 /c 2 — /c 2 , 
which gives 



kl(w) 



uj 2 (uj-N)(uj + N*) 
c 2 (uj - P)(uj + P*) 



^2 _ ,2 
c 2 ' 



(17) 



6 




0.6 0.64 0.68 0.72 

uj/(k x c) 



FIG. 6. The real and imaginary parts of fcaz(w) = 
A v /e20j 2 /c 2 — fc§, with t2 given by Eq. (|16|l . We have set U2 = 
k x c, N = fc x c(6/10 - i/1000) and P = k x c(7/W - i/1000). 
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FIG. 7. The reflected electric field (solid line, z < 0) and 
transmitted electric field (dashed line, z > 0) for a plane wave 
incident to the gain medium described by Eq. (|16[1 and Fig. [6] 
The field is plotted for x — and fe^cf = 10 5 , and normalized 
to the incident field. The amplitude of the reflected field is 
0.98. The parameters used: loi — 0.65k x c, 7 = 10~ 6 , and 
ei = 4. 



Choosing u 2 = k x c, we can tailor the frequency depen- 
dence of k 2z by carefully selecting the locations of zeros 
and poles. Let N = n — iC and P = p — iC, where C > 0. 
All poles and zeros are now located in the (closed) lower 
half- plane. For u> > 0, assuming C « the longitudi- 
nal wavenumber can be written 



k 2 2z {w) = A{ui) (B(Lo)+iC(n-p)) 



(18) 



for real functions A(u>) and B(u>); in addition A(u>) > 0. 
Hence, for n < p, lmk 2z < for all positive frequencies. 
Since k2 Z is analytic in the upper half- plane of u>, and 
since k 2z +^/c as us — > 00, k 2z will be located in the 
forth quadrant of the complex fc 2 z-plane, i.e., Kek 2z > 
and Imfc2z < for all to > 0. A proper evanescent 
or "anti-evanescent" wave has |Re(fc 2z )/Im(fc 2z )| <C 1, 
so we search for values of uj\ satisfying this require- 
ment. Analyzing Fig. [6j there exists an ui\ where 
|Re(fc 2z )/Im(fc2 2 )| <C 1 for n < u\ < p. We have 
hence found a medium for which k 2z describes an "anti- 
evanescent" wave in a finite frequency range. 

Any realistic incident wave contains a spectrum of 
wavenumbers k x . While there are no zeros of k\ z in the 
upper half-plane for the particular k x considered above, 
this is not the case for all possible k x . Thus, also for 
this medium there are growing waves. The fact that 
the medium has large gain, and the presence of insta- 
bilities, mean that it is very challenging to observe the 
"anti-evanescent" response in practice. In principle, how- 
ever, up to a certain time the amplitude of the instabil- 
ities can be limited by ensuring a narrowbanded spec- 
trum of incident k x s. Formally, if a is the width of 
the incident wave, and £ a {x, z, t) is the resulting electric 
field, lim (J _ i . 00 £ a (x, z, t) tends to the "anti-evanescent" 
response as t — > 00, while lim^oo £ a (x, z,t) = 00 for any 
finite a. 

The permittivity (TT5)) has a double pole at cj = 0. 
While the medium is causal in principle, the medium 



might be easier to realize if the pole is moved slightly 
away from the origin, into the lower half-plane. It turns 
out that this modification does not alter the permittivity 
function signifiantly, in the frequency range of interest. 
Also, if desired, the behavior at u> = 00 can be adjusted 
along the lines described in Ref. 17 1. 



In Fig. [7] we plot the reflected and transmitted field 
for the gain medium (|16[) . calculated with the inverse 
Laplace transform for a sufficiently large t when the tran- 
sients have died out. Only a single k x has been excited. 
The reflection amplitude is 0.98, and the transmitted field 
is an exponentially increasing function of z. While real- 
izable in principle, the example is highly unrealistic: To 
observe a behavior similar to that in Fig. t must be at 
least of the order of 10 2 (fca;c) _1 ; otherwise the transients 
would disturb the picture. Any realistic gain medium 
has finite thickness. However, to act as a semi-infinite 
medium, the thickness d of the gain medium must satisfy 
d > ct, or k x d> 10 2 , such that the light has not reached 
the back end. With the "anti-evanescent" growth rate in 
Fig. [71 this would imply unphysically large fields (or in 
practice, nonlinear gain saturation). Hence, if the "anti- 
evanescent" behavior is to be observed experimentally, 
one would need to construct a medium where the tran- 
sients die out rapidly, and/or a medium which leads to a 
sufficiently small | Im k 2z | • At the same time the medium 
must violate the conditions in Sec. that is, it must 
have large gain and/or large dispersion for some frequen- 
cies. 



V. CONCLUSION 

We have considered the case where light is incident 
from a high-index medium to a low-index medium with 
gain, generalizing the situation with total internal reflec- 
tion. 
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In principle, it is apparent that both solutions {k% z in 
the second and fourth quadrant of the complex plane) 
can be attained with a suitably engineered medium. In 
other words, evanescent gain may or may not be the 
case, dependent on the detailed permittivity function. 
This demonstrates the fact that the sign of k2 Z cannot 
be determined from the electromagnetic parameters at a 
single frequency, but must be identified from the entire 
frequency domain dependence, after a check of possible 
non-analytic points (instabilities) in the upper half-plane 
of complex frequency. 

For conventional, weak gain media, we have seen that 
there is an absolute instability associated with infinite 
transverse dimensions. In some cases this instability can 
be eliminated or ignored; then evanescent gain prevails. 



Appendix A: Finite incident beam and finite size 
medium 



To origin from a realistic source, an incident beam 
should not only be causal, but also be of finite width. 
We will here describe how to model an incident beam 
using standard Fourier optics, and argue that even for 
active media, we are allowed to interchange the order 
of integration with respect to transverse wavenumber k x 
and frequency u. Thus we can treat a causal excitation 
of each k x separately. 

Let £(x, t) be the incident TE field at the interface be- 
tween the high-index medium and the active low index 
medium. Performing a Laplace transform t — > uj followed 
by a Fourier transform x — > k x , we obtain the trans- 
formed field E(kx,u). The inverse transform is given by 



£{x,t) 



1 



(2tt 



00+27 



dcue 



(27T) 



— 00+27 
00 pOC 



dk x E(k x ,u)e 



tk j; x 



00 J —00 



E{k x , J + i 7 )e^ i+7 ¥ fc * x dfc x dw'- 

(Al) 



By Fubini's theorem we may interchange the order of 
integration in Eq. (|A1|) . provided E{k x ,oo' + £7) is abso- 
lute integrable with respect to k x and u/. This is the case 
assuming that the incident field is sufficiently smooth 
with respect to t and x. For example, taking the incident 
wave to be a(x)e " lX b(t)e~' lu!xt , the transformed field be- 
comes E(k x , ui' + ij) = A(k x — K x )B(uj' — uji + i'j), where 
A is the Fourier transform of a, and B is the Laplace 
transform of b. Here we assume that b(t) = for t < 0. 
If a and b are continuous, A and B are absolute inte- 
grable. 

We can repeat the above argument for the total field 
(incident + reflected, and transmitted). Assuming no 
superexponential instabilities, the total field is uniformly 
bounded: 




\£{x,z,t)\ < Cexp( 7 i), 



(A2) 



FIG. 8. The semi-infinite gain medium can be replaced by 
a finite size gain medium, provided we only consider times 
t < d/c. 



for positive constants C and 7. Then the transforms 
t — > to followed by x — > k x exist, and we can express the 
total field in the form (|A1|) . The total field is determined 
using the wave equation. In order to consider each mode 
k x separately, we interchange the order of integration for 
each term in the wave equation. To do so, we require the 
second order derivatives with respect to t and x to be 
continuous. 

It remains to prove that our solution is consistent with 
this requirement. From the theory in Sec. [TTl we find the 
solution for each k x , given a sufficiently smooth incident 
field. For this solution, the Fresnel equations show that 
the reflection and transmission coefficients tend to zero 
and unity, respectively, as \uj'\ — > 00 or \k x \ — » 00. There- 
fore the reflected and transmitted field in the (u>,k x )- 
domain adopt any absolute integrability property from 
the incident field. 

In our analysis the incident field u(t) exp(ik x x — iuj\i) 
is not continuous. Hence, strictly speaking, the above 
described method cannot be used. However, by smooth- 
ing the discontinuity around t = 0, we can make the field 
and its second order derivative continuous. This mod- 
ification will not affect the discussion in general, since 
a slower transient will reduce the bandwidth. Thus the 
side waves are excited weaker, such that inequality (|15[) 
is satisfied with an even larger margin. 

In a real experiment, not only the beam width, but 
also the size of the active medium itself, must be finite. 
Provided the fields never reach the end of the structure 
in the time window of interest, the fields will be identical 
to those in a semi-infinite active medium. Thus we can 
consider a setup as in Fig. [8j where the least distance 
from the incident beam to the boundary is d. For t < d/c 
the fields will be the same as if the finite-size medium 
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were replaced by a semi-infinite medium. 



Appendix B: Instabilities in infinite media 

It is convenient to divide instabilities into two cat- 
egories, convective and absolute instabilities (see e.g. 
[20l UM)- Media with absolute instabilities are often 
regarded as impractical for small-signal, linear applica- 
tions, since for an unbounded medium the fields diverge 
even at a fixed point in space. In contrast, media with 
convective instabilities are useful in the linear regime. 
Here the fields do not diverge at a fixed point in space; 
the growing wave is rather convected away. 

However, even in the case with only convective insta- 
bilities, there may be fundamental problems in the case 
where the medium occupies an infinite region or half- 
space: Any small perturbation may propagate an infinite 
distance, thus picking up an infinite amount of gain. In 
our analysis we assume that the active medium is dark 
for t < 0. It is not clear whether this is possible, not 
even in principle, since perturbations in the remote past 
would not die out but rather increase exponentially. 

The remedy is motivated by practical considerations. 
In an experiment, the active medium must have finite 
size in all directions. For a medium without absolute in- 
stabilities and with a given maximum size ci, there will 
be no instabilities provided the gain is sufficiently weak. 
Examples of such configurations include optical ampli- 
fiers, and laser resonators with pumping below threshold. 
When there are no instabilities, we can turn on the pump 
in remote past such that the perturbations have died out 
before t = 0. For < t < d/c we can still regard the 
medium as semi-infinite, since, as seen from Fig. |8l it 
makes no difference. 



Appendix C: Determining the reflected time-domain 

field 



+oo + V), is the same as integrating around all branch 
cuts and poles. The denominator k% z + k 2z does not have 
any zeros, provided the permittivity ei can be considered 
constant and larger than unity in the frequency range of 
interest. Thus we only need to consider the branch cuts 
extending from branch points of ki z and k 2z , and the 
pole at to = u>±. Note that the branch cuts are arbitrary, 
as long as they extend from the branch points. We let 
all branch cuts lie parallel to the imaginary axis, towards 
Imw = — oo. Sec illustration in Fig. [5J The branch 
points of ki z are located far away from (and below) the 
real frequency axis, provided the medium's bandwidth is 
sufficiently large. The wavenumber k 2z has two branch 
points in the upper half-plane, located immediately above 
lu = ±k x c. In addition there are four branch points lo- 
cated below the real frequency axis, with imaginary parts 
— r/2; two simple zeros and two simple poles. The inte- 
grals around the latter four branch cuts decay with time 
constant at most 2/T. Thus, for t > 2/F, the only con- 
tributing terms are the residue of the pole at u)\, and 
the contribution £bc(0,£) from the two remaining branch 
cuts of k 2z : 



£ P (o,t) 



ku 



ku 



■eocp(-iwit) + £bc(0,t). (C3) 



Here k\ z and k 2z have been evaluated at the frequency 
ui\. We write £bc(0,i) = £bc-(0,i) + £bc+(0, t), where 
£bc-(0>£) an d £bc+(0,i) are the contributions from the 
branch cuts in the left and right half-planes, respectively. 

Assuming F <C 1, T < wo and v2wo < k x c, the branch 
cut in the right half-plane extends from approximately 
lu = k x c + iS to to = k x c — ioo, where 8 < FT. Then, for 

* £ 2 / r 



£bc+(o,t) 



l 

1 

27 



p\ (CJ) — : — aw 

k x c i ^ 

Pr(ljj) aw. 



VjJ\ — ILU 

exjp(-iujt) 



k x c—i 



IL0\ — ILU 



(CM) 



Here we will calculate the reflected field in the time- 
domain, when the gain medium is described by a weak, 
inverted Lorentzian function: 



e 2 (w) = 1 - 



iioT 



(CI) 



In Eq. (jCll) F, loq, and T are positive parameters, de- 
scribing the resonance strength, frequency, and band- 
width, respectively. The physical, time-domain reflected 
field at z — x — is given by the inverse Laplace trans- 
form (|13j) . repeated for convenience here: 

E pM = lJ +00+il ^-^l (C2) 



2?r 



— oo+z7 



k\ z + k 2z iuii - iu; 



The field can be interpreted by evaluating integral (IC2[) 
by a generalized version of the residue theorem. We 
here recognize that integrating along path — oo + 17 to 



Here subscripts 1 and r indicate that p(u) is discontinuous 
when crossing the branch cut, denoting the left and right 
side of the branch cut respectively. We further define 
/i,r(w) = k 2z jk\ z . Since k 2z is small in the vicinity of 
k x c, by first order approximation pi jT (u)) = 1 — 2/i jr (u;), 
where / r (w) = — The integral (|C4I) can now be 

simplified: 



£bc+(0,i) 



-2 
in 



^M-pH-ii*,, (C5 ) 



In order to obtain a manageable expression for fi(co), it is 
useful to express k 2z as a function of its zeros and poles. 
With poles denoted by subscript p, and zeros denoted by 
subscripts k x c and loq (indicating the location along the 
real frequency axis), k 2z appears as 



,2 _ ( w ~ ^ ) ( u + ^0 ) ( u ~ u k*c) + <4 aC ) 

/Co- — 



c 2 (lu - cj p )(cj + UJ*) 



(C6) 
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x Pole - Wj 
* Poles - fez 
O Zeros - fo z 



-fee 



bc- 



-COp 



\mco 



cap 



Ci 



k x c 



OJj 



Re® 



-bc+ 



FIG. 9. Circles and stars mark zeros and poles of k-2z respectively. The cross marks the pole at ui = cji. Branch cuts are 
arbitrarily chosen to lie parallel to the imaginary axis, extending into the lower half-plane Im(cj) < 0. The integration path C 
is shown with the dashed line. Contributing branch cuts and poles are enclosed by paths Ci, be— and bc+. For t > 2/T, we 
have § c /(w)dw « § Ci /(w)dw + J bc _ /(w)dw + / bc+ /(w)dw. 



Identifying 5 = Im(wit lC ) and Wj = Im(w), and recogniz 
ing that (u> — uj Uo )/(uj — ui p ) « 1 at u = ^k m c, Eq. (|C6| 
can be simplified: k\ z rj —i2k x {5 — u>i)/c. This gives 



\f } (k x c + iui)\ « y/2(S - WiJ/^cfei - 1)). (C7) 

For fcjC - wi 3> r, we can now find an upper bound of 
integral (|C5]l by noting that - uj 1 < \J FT + Y/2 for 



all Ui considered. We can estimate £bc-(0, t) similarly, 
yielding the bound 



\£bc(0,t)\ < 



T 3/2 



(k x c- u 1 )^k x c(e 1 - 1) 



«c(ei - 1) (' 



e Frt _ e _rt/a 



(C8) 

Consequently for 2/T < t < 1 /FT, the field is well de- 
scribed by the first term in Eq. (|C3|) . 
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